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Introduction 

A second-order unstructured-grid code, developed and used primarily for steady aerodynamic simulations, is 
applied to the synthetic jet in a cross flow. The code, FUN3D, is a vertex-centered finite-volume method orig- 
inally developed by Anderson[l, 2], and is currently supported by members of the Fast Adaptive Aerospace 
Tools team at NASA Langley. Used primarily for design[3] and analysis[4] of steady aerodynamic configu- 
rations, FUN3D incorporates a discrete adjoint capability, and supports parallel computations using MPI. 


Solution Methodology 

A detailed description of the FUN3D code can be found in the references given above. The code is under 
continuous development and contains a variety of flux splitting algorithms for the inviscid terms, two meth- 
ods for computing gradients, several turbulence models, and several solution methodologies; all in varying 
states of development. Only the most robust and reliable components, based on experiences with steady 
aerodynamic simulations, were employed in this work. 

As applied in this work, FUN3D solves the Reynolds averaged Navier-Stokes equations using the one 
equation turbulence model of Spalart and Allmaras[5]. The spatial discretization is formed on unstructured 
meshes using a vertex-centered approach. The inviscid terms arc evaluated by a flux-difference splitting 
formulation using least-squares reconstruction and Roe-type approximate Riemann fluxes. Green-Gauss 
gradient evaluations arc used for viscous and turbulence modeling terms. 

The discrete spatial operator is combined with a backward time operator which is then solved iteratively 
using point or line Gauss-Seidel and local time stepping in a pseudo time. For steady flows, the physical 
time step is set to infinity and the pseudo time step is ramped up with the iteration count. A second-order 
backward in time operator is used for time accurate flows with 20 to 50 steps in the pseudo time applied at 
each physical time step. 

For this effort, FUN3D was modified to support spatially varying boundary and initial conditions, and 
unsteady boundary conditions. Also, a specialized in/out flow boundary condition was implemented to 
model the action of the diaphragm. This boundary condition is described below in more detail. 

The grids were generated using the internally developed codes GridEX[6] for meshing the surfaces and 
inviscid regions of the domain, and for CAD access; and MesherX[7] for meshing the viscous regions. Grid 
spacing in on the surfaces and in the inviscid regions arc indirectly controlled by specifying sources. The 
viscous layers are generated using an advancing layer technique. MeshersX allows the user to control the 
spatial variation of the first step off the surface, growth rates, and the termination criterion by providing 
small problem dependent subroutines. 


Modeling of Diaphragm Boundary 

A specialized in/out flow boundary condition is formulated to model the action of the diaphragm at a station- 
ary boundary. This boundary condition is similar to a standard characteristic boundary condition, commonly 
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applied at far- field boundaries, in that the inviscid boundary flux is evaluated from an intermediate solution 
state I ; that is computed from the current boundary solution I 1, and a prescribed external state U p . In the 
weak form implemented in FUN3D, the boundary solution is not replaced by the intermediate solution but is 
allowed to evolve. The viscous gradients are evaluated using only the boundary and interior solution values, 
and have no knowledge of the prescribed or intermediate solution states. 

At a far-ficld boundary, the intermediate solution is determined by the characteristic relationships that 
model the convective and acoustic waves that arc assumed to exist between the boundary and the far-held. 
For example, at a subsonic boundary 


R+(0) = fl + (0), <r<m = Jr(tW, and JW.) = { (1) 

where h >+ -~' 1 denotes the characteristic variable associated with acoustic waves leaving(entering) the do- 
main, R c denotes characteristic variables associated with convective waves, and V n is the normal velocity 
through the boundary, with outflow taken as positive. 

In the case of a moving diaphragm, the “prescribed” state is that at the diaphragm face, and only the 
diaphragm velocity is known. If the boundary were moving with the diaphragm, than the intermediate state 
would be given by 


v; = 0, R+(Ui) = R+(U b ), and R c (U t ) = R c (U b ) . (2) 

Simply applying these conditions at a stationary location, by setting V n equal to the diaphragm velocity, 
would result in under specifying the flow during the inflow phase of the simulation. To stabilize the simu- 
lation, the latest value of 11 ,, is saved during the outflow phase of the calculation, and reapplied during the 
inflow phase to give: 

V n = ocas(uf), R+(Ui) = R+(U b ), and R c (U t ) = { (3) 

where II, denotes the value saved from the most recent outflow cycle. In practice, only the entropy is saved 
and reapplied. The boundary tangential velocity is allowed to evolve without constraint. Also, during the 
inflow cycle, the entropy is gradually relaxed back toward its initial value. 


Implementation and Case Specific Details 

The geometry of the cavity is simplified by making the diaphragm flat with the height of the diaphragm 
chosen so that the “at rest” volume of the cavity is unchanged. To reduce the problem size, the domain is 
divided at the tunnel centerline and only half of the domain is grided. Although all simulations presented 
here were performed on this half-domain geometry, it is possible to reflect the grid about the centerline to 
obtain a symmetric grid for the larger problem. The computational domain extends from 8 jet diameters 
upstream and to the side of the jet center, and to 16 jet diameters downstream of the jet center. 

The fine grid was sized so that the first spacing on the tunnel wall would have a y+ < 1. Spacings 
inside the cavity were based on preliminary simulations with steady blowing applied at the diaphragm face. 
Griding sources were placed around the lip of the jet and along the jet centerline in an effort to improve 
the clustering there. Figures l(a-c) shows three views of the mesh on the symmetry plane that illustrate the 
mesh distribution near the jet. 

A coarse grid was generated simply by doubling all grid sources and modifying the growth rate of 
the viscous layers. Originally, this grid was intended only to provide rapid turn-around while sorting out 
boundary and initial conditions. The fine grid has 1457853 tetrahedra and 255426 nodes. The coarse grid 
has 254046 tetrahedra and 46063 nodes. 

Characteristic boundary conditions are applied weakly at the inflow, outflow and top boundaries. Sym- 
metry conditions arc imposed at both y=constant boundaries. No-slip conditions arc enforced on the tunnel 
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Figure 1 : Grid on the symmetry plane: a) overall view, b) near jet exit, c) near corner of jet exit. 


floor, the top wall of the cavity, and in the jet contraction. At these no-slip boundaries, the temperature is 
set to the adiabatic wall temperature, and the density is determined from the continuity equation. The side 
walls of the cavity arc treated as inviscid walls, which FUN3D enforces weakly. The in/out flow boundary 
condition described in the previous section is applied on the lower wall of the cavity. 

The inflow and initial conditions in the tunnel region were generated by performing a separate boundary 
layer simulation with FUN3D and extracting the solution at the appropriate Reynolds number. Initially, there 
is no flow in the cavity and the pressure is set to the freestream value. The simulation was stalled impulsively 
with the diaphragm velocity (normalized by u,^) specified as —0.007 cos (u>t). The simulations used 720 
time steps per period, with results saved every 5 degrees of phase. The long time averages requested by the 
workshop were computed by averaging these 5-degree samples. The fine grid simulations required 35 hours 
per period when using 16 intel processors; the coarse grid simulations required less than 1 1 hours per period 
when using 8 processors. 

The fine and coarse grid simulations produced noticeably different solutions in several respects. Fig- 
ures 2 (a) and (b) show the fine grid solution history at a location 0.17mm upstream of the jet center. 
Although the vertical velocity component W settles quickly to a periodic solution, the streamwise velocity 
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Figure 2: Solution at 0.17 mm upstream of jet center: a) and b) fine grid startup history of W and U, 
respectively; c) comparison of U. 


component U requires 4 periods. In the coarse grid simulation, the solution requires even longer, 6 periods, 
before the U velocity component becames periodic. The fine and coarse grid simulations give similar results 
for W, but the U velocity components, compared in fig. 2(c), have completely different character during 
the blowing phase. The coarse grid produces results similar to the experiment, while the fine grid predicts 
a large negative streamwise velocity component during the blowing cycle. The contour plots of U, shown 
in fig. 3, indicate that this negative streamwise velocity is due to a vortical behavior that develops in the jet 
exit flow. 
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Figure 3: Streamwise velocity contours from fine grid simulation during peek blowing: phase = 160. 
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